Circularization and Final Spin in Eccentric Binary Black Hole Inspirals 
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We present results from numerical relativity simulations of equal mass, non-spinning binary black 
hole inspirals and mergers with initial eccentricities e < 0.8 and coordinate separations D > 12 M of 
up to 9 orbits (18 gravitational wave cycles). We extract the mass Mi and spin at of the final black 
hole and find, for eccentricities e < 0.4, that Of/Mf « 0.69 and Mf/Madm ~ 0.96 are independent of 
the initial eccentricity, suggesting that the binary has circularized by the merger time. For e ^ 0.5, 
the black holes plunge rather than orbit, and we obtain a maximum spin parameter at /Mi « 0.72 
around e = 0.5. 
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The field of numerical relativity (NR) has now en- 
tered a stage where binary black hole (BBH) simulations 
can reliably be used to investigate a vast range of inter- 
esting phenomena. Studies have produced gravitational 
waveforms from binary systems in essentially circular or- 
bits [E 0) H, I3|> involving spinning black holes as well 
as unequal mass systems. The level of numerical accu- 
racy achieved by these codes is impressive [3, Q , and in 
some of these studies, the initial binary separations were 
such that it was feasible to directly compare with post- 
Newtonian (FN) waveforms 0, 0]. Other examples of 
exciting new results in NR are investigations of the kick 
imparted to the final black hole (BH) 3, 'e', II, i, 0, [H , 
the spin dynamics of the merging BHs 12, 13] , and, of rel- 
evance to the present work, the m erg er threshold between 
bound and unbound BBHs [3, [l^, and studies of the 
result of a circular inspiral of spinning BHs [H, [13, EBl ■ 
All of this has been possible since the pioneering work 
of Briigmann et al. fl9| . Fretorius [20|, Baker et al. 
[2]| . Campanelli et al. [22|. 

It is well known that gravitational radiation leads to 
circularization of a binary system [23j . In this work, we 
study this circularization in the nonlinear regime. In 
Ref. [T'], it was found that, for equal-mass, non-spinning 
BHs initially in quasi-circular orbits, the merger pro- 
duced a BH with spin parameter af/Aff 0.69, which, 
within the accuracy of the results, was independent of 
the initial separation. 

In this Letter, our main goals are (1) to investigate 
whether sufficient eccentricity is lost during the late 
stages of inspiral to circularize the orbit and exhibit the 
same universality as in the circular case and (2) to ex- 
tract the spin parameter and mass of the final BH, and 
compare the values with those from circular inspirals. 

Although isolated stellar mass BBHs will have com- 
pletely circularized by the time they are observable by 
ground-based interferometers, scenarios have been sug- 
gested for which BBHs in eccentric orbits are not only 
astrophysically interesting but also could be detected by 
space- or ground-based interferometers [l^, [1^ . For in- 
stance, galactic mergers leave behind massive BBHs that 
likely interact with a gaseous environment. A gaseous- 
gravitational driven inspiral could yield a BBH arriving 
at the last few orbits and merger with a non-vanishing 
eccentricity. An observation of the gravitational waves 
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FIG. 1: Waveform polarizations /i+ for the cases e = 0.1 and 
e = 0.3 



from an eccentric BBH merger will allow us to determine 
the amount of angular momentum lost to gas and, in 
particular, the gravitational torques between the binary 
and a circumbinary disc that affect the eccentricity of the 
binary [26!|. 

Methods: We construct initial data using the puncture 
approach [23|, which requires specifying the coordinate 
locations and momenta for the two BHs. For "circular" 
orbits, we follow HI]. For eccentric models, we use the 
conservative 3FN expressions in Ref. [1^. These expres- 
sions require the specification of the eccentricity e and 
the mean motion n = 21: / Pr, where Pr is the radial (peri- 
center to pericenter) orbital period. There are three FN 
eccentricities, which are the same to 1 FN order, and we 
choose Cf, which appears in the FN Kepler equation, fol- 
lowing Ref. [29| . It is important to keep in mind that the 
eccentricities we quote (and we use them also to label the 
models) are to be taken only as a guide to the eccentric- 
ity in the initial data, as the FN expressions used do not 
include radiation reaction, and the FN parameters are 
in a different coordinate system to the puncture initial 
data. 

We construct a family of initial data by fixing n — 
0.01625/M {Pr 387M) and varying e in the range 
0.05 — 0.8 (note that to 2FN order, this means that the 
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TABLE I: Initial data parameters: The runs are labeled by 
their initial eccentricity e. The BHs have hnear momenta 
±Pi,2/M and are separated by a coordinate distance D/M. 



systems have the same binding energy and that, at high 
eccentricities, there are portions of the orbit for which the 
PN condition ii/c <C 1 is no longer valid). The bi nary 
separation D is determined from Eq. (23) in Ref. [29j . 
and the tangential linear momentum, P/M, of each BH 
at apocenter is obtained from J = PD, where J is the 
total angular momentum computed as a PN expansion in 
n and e (Eq. (21) in Ref. [29]). The bare BH masses mi^2 
are chosen to make the irreducible BH masses Mi^2 = 0.5 
(i.e. M = Ml + M2 = 1). Table H] provides the initial 
data parameters. 

The numerical simulations and results in this work 
were obtained with the same infrastructure used in our 
previous BBH studies (see Ref. fs^l for full details). We 
have evolved the circular model at three different res- 
olutions (finest grid spacings of M/38.7, M/51.6 and 
M/64.5). We obtain approximately fourth order con- 
vergence in the total energy and angular momentum ra- 
diated, consistent with the designed 4th order accuracy. 

Results: In Figs. [T] and [5] we display the gravitational 
wave strains and coordinate inspiral tracks for e = 0.1 
and e = 0.3. It is evident that the difference in ini- 
tial eccentricity has a large effect during the inspiral. 
Qualitatively, the case with larger eccentricity exhibits 
a more rapid inspiral ^23*1. However, at some point both 
systems enter a "circular" plunge, hinting that circular- 
ization may have occurred. We find that the simulations 
with e > 0.5 show plunge- type rather than orbital- type 
behavior in the coordinate motion from the very start. 
Note that the tracks shown in Fig. [5] represent the coordi- 
nate positions of the individual BHs, and once a common 
horizon forms, they are less meaningful. 

We now consider the emitted radiation and focus on 
the dominant £ = 2, m = 2 mode of the complex 
Newman-Penrose (NP) quantity \E'4 = A{t) exp {—iip{t)). 
To compare the orbits, we apply a time shift to A and 
(/?, so that the maximum of A (i.e. the peak of the am- 
plitude of the gravitational wave) is at t/Mi = in each 
simulation. In Fig. [31 we plot the shifted amplitudes and 
frequencies u; — d(p/dt extracted at r = TOM. The cases 
displayed are those with eccentricities e = — 0.5 in steps 
of 0.1 and e = 0.8. 

In Fig. [31 the oscillations and growth in uj at early 
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FIG. 2: Inspiral tracks for initial eccentricity e = 0.1 (left 
panel) and e = 0.3 (right panel). 



times (inspiral) can be in general terms understood from 
simple Newtonian considerations. That is, ignoring radi- 
ation reaction, the oscillations (i.e. amplitude and period) 
in uj are a direct consequence of the eccentricity and not 
present in the e = case. The period of these oscillations 
is the period P,. of the extrema in the separation, and the 
amplitude of the oscillations increases with e. The addi- 
tion of radiation reaction leads to an overall growth of lo 
with time due to the energy and angular momentum loss, 
and this is clearly visible in the figure. The amplitude 
of the oscillations in co should decrease with time, cor- 
responding to a reduction in eccentricity. However, over 
these timescales, it is difficult to separate this effect from 
the secular increase in lo. Also consistent with the pre- 
dictions in [23], the higher eccentricity evolutions merge 
more quickly. We note that in the Newtonian case, we 
would observe that Pr — P^, where P^ is the time the 
binary takes to complete one revolution in the angular 
coordinate (j). Due to the effects of precession caused by 
general relativity, the two periods are very different (this 
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FIG. 3: Frequency of the ^ = 2, m = 2 mode of the NP 
radiation scalar r ^'4 
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can be seen from PN equations). 

During the merger or plunge phase, lo increases dra- 
matically and then levels off, signaling that the BBH has 
merged. After this point, tu remains constant, a direct 
consequence of the quasi-normal mode (QNM) ringing of 
the final BH. 

As mentioned before, one of the objectives of this work 
is to investigate whether a given initial data configura- 
tion will circularize before it merges. By this we mean 
that the radiation from the late stages of the evolution 
is identical to that from an orbit which started with zero 
eccentricity. We see in Fig.[3]that the low e evolutions ap- 
proach the same final state as a circular orbit at the very 
late stage of inspiral. The lus from different eccentrici- 
ties near t/Mf = seem to be indistinguishable for low 
enough eccentricity. In order to investigate this in more 
detail, in the inset of Fig. [3] we focus on the plunge stage. 
Here we plot eccentricities e = — 0.8 in steps of 0.1. Up 
until e = 0.4 and after t/M[ « —50, the frequencies uj 
from each run follow each other. Noticeable differences 
start showing for e > 0.5, which is the first configuration 
to plunge immediately without orbiting first. 

We now discuss Aff and at, computed using two in- 
dependent methods. In one method, they are obtained 
from the radiated energy and angular momentum using 

Aff = Madm - -Ei-ad and af/M{ = (Jadm - Jrad)/Mf. In 

the second method, Mf and af are computed from the 
QNM frequencies [31| emitted by the final BH, extracted 
using least squares fitting. As a cross-check, for some of 
the models we also determine Of/Mf using an approxi- 
mate technique derived from the isolated horizon formal- 
ism [l^ [s^] • Table |TT] gives the energy i?i-ad and angular 
momentum Jjad radiated as well as the final mass Mf 
and spin Of. Figure S] gives the final mass Mf and spin 
Of as a function of e. Notice the agreement in af and Mf 
that the three methods give within the estimated error 
bars. The final mass and spin also agree well in the circu- 
lar case with the values obtained in Ref. [s^l- The error 
bars for the QNM-derived quantities are dominated by 
the uncertainties in the fitting procedure, which are esti- 
mated as the variations of the fit parameters over a range 
of fitting windows. The errors on the radiation-balance 
quantities are dominated by the finite differencing error 
of the simulations. Due to excessive computational ex- 
pense, we have not run very high resolution versions of 
the eccentric simulations, and so use the errors from the 
corresponding low resolution circular orbit as a rough 
guide to the errors in the eccentric cases. 

Given an initial eccentricity, it is possible to choose a 
large enough semimajor axis or orbital period for which 
the binary circularizes before it arrives at the merger. 
Our family of initial configurations was designed to in- 
vestigate, for a fixed initial orbital period, how much ini- 
tial eccentricity a binary is able to have and still enter 
the merger with essentially vanishing eccentricity. Since 
we do not have a good measure of eccentricity appli- 
cable prior to the merger, we focus on the end state, 
namely Mf and af of the final BH. We see from Fig. 0] 
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TABLE II: Extracted quantities: Energy i?rad and angular 
momentum Jrad radiated; final spin parameter at and mass 
Mf computed from J^ad and -Erad as well as from QNM ring- 
ing. 
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FIG. 4: Plots of Mf/Madm and at/Mi as functions of the 
initial eccentricity e. Note that the error bars shown here 
on the radiation quantities for eccentric runs are taken from 
the low resolution circular case and thus should be treated as 
indicative only. 



that af/Mf « 0.69 for e < 0.4 and Mf/Madm ~ 0.96 for 
e < 0.5, both values of A/f and af in agreement with the 
circular result. We note that the remaining orbits, the 
ones which do not circularize, are all configurations which 
seem to plunge immediately rather than entering an or- 
bital phase. We conclude that for the systems we studied 
with approximately constant initial orbital period, within 
our error bars, orbits with e < 0.4 essentially circularize 
before they merge, and orbits with e ^ 0.5 plunge. 

We also observe for e ^ 0.4 that rather than af decreas- 
ing monotonically, a maximum spin parameter af/Mf « 
0.72 is obtained around e = 0.5. Given the size of our 
uncertainties and that the maximum is found in the three 
independent methods used to calculate the spin, we are 



4 



confident that this maximum is real for our family of ini- 
tial data. At about e — 0.6, Of starts decreasing mono- 
tonically. We are currently considering larger, but still 
computationally feasible, initial separations to investi- 
gate if there is any bound orbital (rather than plunge) 
configuration that docs not circularize. 

As e — > 1, corresponding to vanishing linear mo- 
menta (i.e. a head-on collision from rest), we find that 
at/Mf — !■ 0, in line with the symmetry of the head-on 
collision, and Mf ~ Madmj as expected, since NR sim- 
ulations of a head-on collision have shown that Mf ^ 
(1 - 0.001)Madm [13. Note that the ringdown result for 
e = 0.8 gives Mf > Madm which is clearly unphysical, 
but the error bars account for this. 

Conclusions: We have carried out a series of eccen- 
tric orbit simulations of BBH systems in full nonlinear 
general relativity to investigate the merger regime and 
final BH. The family of simulations consisted of bina- 
ries with approximately constant initial orbital period 



and varying initial eccentricity. We found that for initial 
e < 0.4, the final BH parameters are Aff/A/adm ~ 0.96 
and af/M[ « 0.69, the same as in the circular case. As 
a consequence of this, we also found that for e < 0.4 the 
binary begins to enter a universal plunge at t 50 Mf be- 
fore the amplitude of the gravitational radiation reaches 
its peak. 

While preparing the manuscript of this work, a study 
by Sperhake et al. J35| appeared with both similar and 
complementary conclusions to those in our present work. 
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